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Abstract 

t> ' Aiming to observe cosmic gamma rays in the 10—1000 TeV energy region, we 

propose a 10000 m 2 underground water-Cherenkov muon-detector (MD) array that 
operates in conjunction with the Tibet air-shower (AS) array. Significant improve- 
ment is expected in the sensitivity of the Tibet AS array towards celestial gamma- 
■ ray signals above 10 TeV by utilizing the fact that gamma-ray-induced air showers 

contain far fewer muons compared with cosmic-ray-induced ones. We carried out 
detailed Monte Carlo simulations to assess the attainable sensitivity of the Tibet 
AS+MD array towards celestial TeV gamma-ray signals. Based on the simulation 
results, the Tibet AS+MD array will be able to reject 99.99% of background events 
at 100 TeV, with 83% of gamma-ray events remaining. The sensitivity of the Ti- 
bet AS+MD array will be ~ 20 times better than that of the present Tibet AS 
array around 20—100 TeV. The Tibet AS+MD array will measure the directions 
of the celestial TeV gamma-ray sources and the cutoffs of their energy spectra. 
Furthermore, the Tibet AS+MD array, along with imaging atmospheric Cherenkov 
telescopes as well as the Fermi Gamma-ray Space Telescope and X-ray satellites 
such as Suzaku and MAXI, will make multiwavelength observations and conduct 
morphological studies on sources in the quest for evidence of the hadronic nature 
of the cosmic-ray acceleration mechanism. 
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1 Introduction 



Supernova remnants (SNRs) are theoretically considered the most plausible 
candidates for acceleration of hadronic cosmic rays up to the knee at ~ 10 15 eV 
in the cosmic-ray energy spectrum. 7r° decays following inelastic collisions 
between those accelerated charged cosmic rays and ambient matter, such as 
molecular clouds near the acceleration sites, are expected to produce gamma 
rays with energies as high as 100 TeV. Meanwhile, accelerated electrons in 
SNRs have difficulty producing gamma rays with energies above 100 TeV via 
bremsstrahlung or inverse Compton scattering, due to synchrotron cooling by 
magnetic field and the reduction of the cross section by the Klein-Nishina 
effect. Observation of the acceleration limit of cosmic gamma rays, therefore, 
is important to establish the acceleration of hadronic cosmic rays in SNRs. 

The H.E.S.S. group reported on the discovery of 14 new TeV gamma-ray 
sources in the southern hemisphere in 2006 ([l|). Their energy spectra were 
measured above 200 GeV and were found to extend up to around 10 TeV, 
many of which have relatively hard spectral indices at TeV energies ranging 
from —1.8 to —2.8. This implies proton acceleration in the sources, since shock 
acceleration models (0; 0) favor ~ —2 as the indices of gamma-ray differential 
energy spectra at TeV energies, although no definite conclusions have been 
drawn yet. Another feature of these sources is their faintness in other wave- 
lengths, which strongly stimulates unbiased celestial surveys independent of 
observations in other wavelengths using wide field-of-view apparatuses, such 
as air-shower arrays. 

Uchiyama et al. reported on the detection of small-scale X-ray flares in a 
supernova remnant RX J1713— 3946, a prominent TeV gamma-ray source in 
the southern sky in 2007 They claimed that these flares suggest an en- 
hancement of magnetic fields in the remnant up to ~ 1 mG, two orders of 
magnitude stronger than normally assumed. For such strong magnetic fields, 
the X-ray emission would be produced by a small population of electrons, 
which in turn could not account for the observed TeV flux of the remnant. 
Thus, they concluded that this provides a strong argument for acceleration of 
protons and nuclei to energies of 1000 TeV and beyond. This statement, how- 
ever, was subsequently rectified by Y. Butt et al. (|5|), who contended that the 
speculated strong magnetic fields and the hadronic origin of the TeV gamma- 
ray emission would lead to radio synchrotron emission ~ 1000 times brighter 
than any previously observed, if the number ratio of electrons to protons is 
not be unrealistically small: < 10~ 6 . They pointed out as well that even if 
~ 1 mG magnetic fields are acceptable, they could not account for the TeV 
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gamma-ray emission present over the entire SNR shell, insofar as they were 
confined within minuscule X-ray hot spots. Tanaka et al. (@) claimed that the 
multi-wavelength spectrum of RX J1713— 3946 based on the Suzaku X-ray and 
the H.E.S.S. TeV gamma-ray data hardly can be explained within the pure 
leptonic model, while the hadronic model with a strong magnetic field provides 
a good fit to the spectrum, with the special arrangements of the model pa- 
rameters being necessary to explain the lack of thermal X-ray emission. Drury 
et al. (7) recently showed that the heating of the ion component of plasma 
caused by the particle acceleration in strong nonlinear shocks can be strongly 
suppressed, so that thermal X-ray emission may no longer be expected. Un- 
der these controversial circumstances, evidence for the hadronic nature of the 
particle acceleration in RX J1713— 3946 remains far from conclusive. 

From this discussion, one can see that observations of the energy spectra of 
TeV gamma-ray sources above 10 TeV, along with morphological studies on 
the sources and positional comparison between their radio, X-ray, GeV, and 
TeV images, are highly encouraged at this moment in order to obtain observa- 
tional proof of the theoretical conjecture for hadronic cosmic-ray acceleration 
in SNRs. An air-shower array with a high background cosmic-ray rejection 
power that observes cosmic gamma rays in the 10—1000 TeV energy region 
at high altitude is considered suitable for this aim. There has been, however, 
no such experiment so far. In this paper, we propose a 10000 m 2 underground 
water- Cherenkov muon-detector (MD) array that operates in conjunction with 
the Tibet air-shower (AS) array. 



The Tibet air-shower array is located at Yangbajing (90.522° E, 30.102° N, 
4300 m above sea level) in Tibet, China (@;S;|3). This array has been suc- 
cessfully operating, leading to a good deal of data on gamma-ray astronomy 
as well as on cosmic-ray physics: multi-TeV gamma rays from acti ye g alac- 
tic nuclei in the flaring state such as Mrk 501 (fiol ) and Mrk 421 (11), the 
steady emission of gamma rays from the Crab Nebula (0; 12), the cosmic-ray 
anisotropy in the sky (1131 ; Il4l ). the cosmic-ray energy spectra around the knee 
(flol ). and so on. In particular, it may be worthwhile to stress here that a few 
spatially separated, small-scale enhancements were observed in the Cygnus 
region at multi-TeV energies (jl6l ). suggesting small-scale (~ 2°) excesses of 
gamma rays due to some extended gamma-ray sources. Actually, following 
this result, MILAGRO reported the detection of TeV gamma-ray sources in 
the Cygnus region (fl7l ). 



Adding the MD array to the Tibet AS array will significantly improve the 
sensitivity towards celestial gamma-ray signals above 10 TeV, by utilizing the 
fact that gamma-ray-induced air showers contain far fewer muons compared 
with cosmic-ray-induced ones. The detailed structure of the MD array and 
the procedure of our simulation study to deduce its attainable sensitivity are 
described. The sensitivities of several present and future experiments are com- 



3 



pared as well. 



2 Tibet Air-Shower Array and Muon Detector Array 



Fig. [T] shows the currently proposed configuration of the Tibet AS array 
with the MD array (the Tibet AS+MD array). The surface AS array (1081 
scintillation detectors, 83000 m 2 ) is a simple extension of the current Tibet air- 
shower array (789 scintillation detectors, 37000 m 2 ). The MD array consists 
of 12 pools, each with 16 cells, covering a total area of 9950 m 2 . Each cell of 
the MD array, shown in Fig. [2j is composed of a concrete water tank 7.2 m 
wide, 7.2 m long, and 2.4 m deep, equipped with two downward-facing 20- 
inch-in-diameter (20" 0) PMTs (HAMAMATSU R3600) on the ceiling. The 
inner surface is coated in white paint so that it acts not only as a waterproof 
layer but also as a reflector of Cherenkov photons, which are emitted by air- 
shower muons and subsequently collected by the PMTs. In order to suppress 
air-shower electromagnetic components, the MD array is set up beneath a 
2.0 m thick layer of soil. The mass thickness of the soil and the concrete 
ceiling (0.5 m) is equivalent to ~ 19 radiation lengths altogether. The energy 
threshold for muons is approximately 1 GeV. 



Our Monte Carlo simulations prove that the configuration of the MD cells 
as described in Fig. [1] provides a sensitivity enhancement of the Tibet AS array 
almost equal to that achieved by an ideal configuration, where the MD cells 
are uniformly distributed in the underground of the AS array; the sensitivity 
difference would be <~ 5%. Considering cost performance and construction 
feasibility, the Fig. [T] configuration has been chosen. The two vertical pil- 
lars reinforcing the cell structure (see Fig. [2]) have only an insignificant effect 
(< 10%) on the amount of the Cherenkov light detected. Water- Cherenkov 
detectors significantly exceed scintillation detectors in cost performance as 
well as in the ability to suppress the influence of environmental background 
radioactivity owing to the energy threshold for the production of Cherenkov 
photons. 



The background noise of the MD array is predominated by relativistic ac- 
cidental muons arriving at a rate of approximately 15 kHz per MD cell (Il8l ). 
This is a much higher rate than 100 Hz, the dark pulse rate of a 20" (f) PMT 
measured in the laboratory at a gain of 10 7 under a supplied voltage of 1500 V 
and at a threshold of 10 photoelectrons. Setting a time window of 200 ns in 
offline analysis reduces the event rate of accidental muons down to 0.6 per air 
shower (15 kHz x 16 cells x 12 pools x 200 ns). 
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3 Monte Carlo simulations 



Air showers are generated with the Corsika Ver.6.204 code (Il9l ). using 
QGSJETOlc as a hadronic interaction model. Primary gamma rays are gen- 
erated along the Crab's orbit with a differential power-law spectrum of E~ 2S 
above 0.3 TeV. Primary cosmic rays, the chemical composition model of which 



(ll5l ) is based mainly on direct observational data, are also generated along the 
Crab's orbit and in the same energy range. The number of generated air show- 
ers is 3.9 x 10 7 (1.8 x 10 8 ) for gamma rays (cosmic rays). 

The generated air showers are thrown at random positions within 300 m 
from the center of the AS array. After simulating the resp onse of the AS 



array for the air showers with the Epics uv8.00 code (1201 ) . we obtain the 
estimated/true arrival direction, the core position, the sum of particle densities 
for all the detectors (Hp), and so on for each air shower. The air showers that 
trigger the AS array whose cores hit more than 300 m from the center of the 
AS array are estimated to be less than 0.8%, 0.04%, and 0.01% of all the air 
showers triggering the AS array at 32 < P < 56, 100 < J2p < 215, and 
215 <J2p< 464, respectively. We select air showers that fulfill the following 
criteria: 

(i) Any four-fold coincidence must occur in the scintillation counters recording 
1.25 particles or more in charge. 

(ii) Five out of the six hottest counters must be contained in the fiducial area 
enclosed by a dashed line in Fig. [TJ For air showers that hit fewer than 6 
counters, all of the hit counters must be in the fiducial area. 

(iii) The residual error in reconstructing the arrival direction by means of the 
least-squares method must be less than 1.0 m. 

In point-source analyses, we normally extract air-shower events coming from 
within a search window around a target source. The radius of the window is 
expressed as Q.9°/y/Hp as a function of Hp, which maximizes the detection 



significance of the sources (11 II). Correspondingly in this simulation analysis, 
we extract events such that the deviation of their estimated arrival direction 
from their true arrival direction is < 6.9°/ \/Hp. 

The trigger efficiency for gamma rays and cosmic rays reaches 100% above 10 
and 50 TeV, respectively. The angular and energy resolutions of air showers 
are estimated to be 0.2° and 40%, respectively, for 100 TeV primary gamma 
rays. 

To simulate the response to the MD pool cells as well as that to the overbur- 
dening, which consists of 2.0 m thick soil and 0.5 m thick concrete, we adopt 
GEANT4 8.0 code (]2 If) , taking into account the details of the cell structure 
shown in Fig. [2j The secondary particles of the remaining air showers are 
traced in the soil, which is assumed to have a mass density of 2.0 g/cm 3 and 
is made up of 70% Si0 2 , 20% A1 2 3 , and 10% CaO. All the particles surviv- 
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ing under the soil layer with energies above their Cherenkov thresholds are 
traced further in the MD array, along with the trajectories of Cherenkov pho- 
tons that they emit. The concrete frame of the MD cells is assumed to have 
a density of 2.3 g/cm 3 and to be made up of the same compositions as the 
soil. The inner surface of the cells is assumed to reflect Cherenkov photons 
isotropically at 70% probability. Fig. [3] shows the quantum efficiency of 20" 
PMTs employed in the simulation. The attenuation length of light in water 
is set to be wavelength-dependent, being 40 m at 400 nm. In the end, the 
numbers of photoelectrons averaged over two PMTs in every MD cell (Ap E ) 
are summed up, taking into account the light collection efficiency of PMTs 
(~60%) including the geomagnetic field effect (1221 ) . The typical value of iVpE 
for vertically penetrating muons that uniformly hit the ceiling of an MD cell 
is estimated to be 17l^%° photoelectrons, as shown in Fig. HI 



4 Results and discussion 



Using the air-shower events that survive the event selection criterion de- 
scribed in §3] and remain within the search window with a radius of 6.9°/v / Ep, 
we obtain the distribution of £iVp E as a function of Ep as shown in Fig. [5j 
EiVp E denotes the sum of iV PE s from the MD pool cells that record photo- 
electrons corresponding to more than 0.15 muons, whereas Ep represents the 
sum of particle densities for all the AS scintillation counters. Ep = 1000 cor- 
responds approximately to 100 TeV for primary gamma rays. EAp E includes 
the noises induced by accidental muons (see $2]) as a Poisson distribution with 
an average of 0.6 events per air shower. At low energies < 10 TeV, more 
than 50% of the gamma-ray events do not hit the MD array (EVp E = 0), 
and among those with EAp E ^ 0, more than 90% are due to contamination 
of the accidental muons. The EAp E -based event-selection criterion, shown by 
the solid line in Fig. [5], is set so as to maximize the detection significance 
of a point source, taking into account the fluxes of the gamma rays and the 
background cosmic rays under the assumption of the one-year observation of 
the Crab Nebula. Selecting air-shower events under the solid line leads to 
the survival efficiencies shown in Fig. [7J Note that the survival efficiencies 
defined here correspond to the event selection by the MD array only, and 
that they do not include the event selection by the search window, which de- 
pends on the angular resolution. The background events are rejected down 
to the 2 x 10- 3 (5 x 10~ 5 ) level at 10 TeV (100 TeV), with 60% (83%) of 
the gamma-ray events remaining. Fig. [S] shows the expected Ep spectra of 
the gamma-ray and background events before and after the event selection by 
EA^p E as shown in Fig. [5j assuming the one- year observation of the Crab Neb- 
ula with a flux of F G = 6.45 x 10~ 10 x (£/0.3TeV)~ 2 - 6 [cm-V^eV" 1 ]. The 
integral flux of the background cosmic rays is assumed to be F CR (> 0.3TeV) = 
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1.03 x 10 4 [cm 2 s x sr x ]. One can see that the number of background events 
above ~50 TeV is well suppressed down to <1. 

The integral flux sensitivity of the Tibet AS+MD array towards a point-like 
gamma-ray source is calculated as shown in Fig. [91 One can see that in one 
calendar year it will achieve 5a sensitivities of 7% and 20% of the Crab flux 
above 20 and 100 TeV, respectively. Note that its sensitivity above 100 TeV, 
where the number of background cosmic rays is fully suppressed down to ~ 1 
event/(l° x 1°), is defined as a flux corresponding to 10 gamma-ray events. 

Table [T] compares the performance of the Tibet AS+MD array with the 
CASA-MIA experiment (|26l ). The angular resolution of the Tibet AS array 
for primary gamma rays (0.2°) is 10 times better than that of CASA-MIA 
at 100 TeV. Therefore, we can take a search window with a 10 times smaller 
radius, so that the solid angle of the window becomes 10 2 times smaller. The 
effective area of the upgraded Tibet AS array (83000 m 2 ) is 1/3 that of CASA- 
MIA, and the background survival efficiency due to event selection by the MD 
array (5 x 10 -5 ; see Fig. [7]) will be 400 times that at 100 TeV. Thus, we can 
make a rough estimation that the sensitivity of the Tibet AS+MD array will be 
y/lO 2 x (1/3) x 400 ~ 100 times better than that of CASA-MIA at 100 TeV. 
This excellent sensitivity may surpass that of existing imaging atmospheric 
Cherenkov telescopes above 12 TeV and the future Cherenkov Telescope Ar- 
ray (CTA) project above 40 TeV as well as the High Altitude Water Cherenkov 
(HAWC) experiment above 15 TeV. The Tibet AS+MD array will be able to 
detect ~10 sources such as the Crab Nebula, TeV J2032+4130, Mrk 421, and 



the MILAGRO sources (lUj) 



In order to demonstrate the sensitivity towards spatially extended sources, 
we take MGRO J2019+37 as an example and analyze the simulation data 
using a search window with a radius of 1.7° independent of Sp (cf. §3]) after 
dispersing the air-shower events randomly within a radius of 0.32° around 
the point source, in accordance with the MILAGRO observation (1171 ). Fig. 
[TU1 shows that the Tibet AS+MD array has advantages of a high signal-to- 
noise ratio and a wide field of view in observing spatially extended gamma-ray 
sources. 



Fig. [TT] shows the sensitivity of the Tibet AS+MD array towards diffuse 
gamma rays from the Galactic plane. The Tibet AS+MD array will be able to 
detect the inverse Compton component from the inner Galactic plane within 
a year. Even if the inverse Compton hypotheses were wrong, the 7r°-decay 
component would be detectable at ~30 TeV within a few years. The inverse 
Compton component from the outer Galactic plane is likely to be detectable 
within a year unless the electron injection spectral index is too soft (~ —2.4). 
Positive detection of the 7r°-decay component from the outer Galactic plane 
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would necessitate several years' observation. Measurements of the flux and 
spatial distribution of these gamma rays in the 10—1000 TeV energy region 
will place constraints on the model predictions, and thus make it possible to 
investigate the relative contribution of the hadronic and leptonic mechanisms. 

Fig. [TJ] shows that the Tibet AS+MD array would be able to detect most 
of the 14 H.E.S.S. sources reported in y]), if located at the H.E.S.S. site. The 
H.E.S.S. experiment has discovered 46 new TeV gamma-ray sources thus far 



(140l) . The Tibet AS+MD array could discover some unknown TeV gamma- 
ray sources in the northern sky as well, as no extensive searches for TeV 
gamma-ray sources have been carried out thus far in the northern sky with 
high-sensitivity apparatuses comparable to H.E.S.S. 

Lastly, multi-TeV gamma rays from intense gamma-ray bursts may be de- 
tectable, taking into account the absorption by intergalactic infrared photons 



(14 ll). if they arise in the relative vicinity of the Earth (z <~ 0.1) and if their 



keV-MeV spectra extend to higher energies. 

The construction of the MD array is estimated to cost roughly 5 million 
USD. The Tibet AS+MD array will be a cosmic-ray experiment designed 
to observe 10—1000 TeV celestial gamma rays and involving the operation 
of the world's largest water- Cherenkov muon-detector array in conjunction 
with the surface air-shower array. This experiment will be complementary to 
H.E.S.S. and the proposed CTA project, both of which are located in the 
southern hemisphere, as well as to MAGIC, VERITAS, and the proposed 
HAWC experiment, which observe up to ~ 10 TeV. The Tibet AS+MD array, 
along with imaging atmospheric Cherenkov telescopes as well as the Fermi 
Gamma-ray Space Telescope (10 keV— 300 GeV) and X-ray satellites such as 
Suzaku (0.4-600 keV) and MAXI (0.5-30 keV), will make multiwavelength 
observations and conduct morphological studies on sources in the quest for 
evidence of the hadronic nature of the cosmic-ray acceleration mechanism. 



Acknowledgments 



We thank the members of the Tibet AS7 collaboration for valuable discus- 
sions. T.K. Sako is supported by a Grant-in-Aid for JSPS Fellows, No. 20-5715. 
For this Monte Carlo calculation we used the computers in the Institute for 
Cosmic Ray Research, the University of Tokyo. 



References 

[1] F. Aharonian et al., Astrophys. J. 636 (2006) 777. 



8 



2 
3 
4 
5 
6 
7 
8 
9 

10] 

11] 
12] 

13] 

14] 

15] 

16] 

17] 

18] 

19] 

20] 

21] 



E. Berezhko, & J. Volk, Astropart. Phys. 7 (1997) 183. 
E. Berezhko, & J. Volk, Astronomy & Astrophysics 357 (2000) 283. 
Y. Uchiyama et al., Nature 449 (2007) 576. 
Y. Butt et al., MNRAS 386 (2008) L20. 
T. Tanaka et al., Astrophys. J. 685 (2008) 988. 
L. Drury et al., Astronomy & Astrophysics 496 (2009) 1. 
M. Amenomori et al., Phys. Rev. Lett. 69 (1992) 2468. 
M. Amenomori et al., Astrophys. J 525 (1999) L93. 
M. Amenomori et al, Astrophys. J 532 (2000) 302. 

Astrophys. J 598 (2003) 242. 
astro-ph/0810.3757vl, accepted by Astrophys. J. 
Phys. Rev. Lett. 93 (2004) 061101. 
Astrophys. J 626 (2005) L29. 
Astrophys. J 678 (2008) 1165. 
Science 314 (2006) 439. 
A. Abdo et al., Astrophys. J 664 (2007) L91. 
M. Amenomori et al., Astrophys. Space Sci. 309 (2007) 435. 
D. Heck et al., Report FZKA, 6019, Forschungszentrum Karlsruhe, 1998. 
K. Kasahara, |http : / / cosmos . n . kanagawa-u .ac.jp/ EPICS Home | 
S. Agostinelli et al., Nucl. Instrum. Methods Phys. Res. A 506 (2003) 



M. Amenomori et al., 
M. Amenomori et al., 
M. Amenomori et al., 
M. Amenomori et al., 
M. Amenomori et al., 
M. Amenomori et al., 



250. 



22] 
23] 
24] 
25] 
26] 
27] 
28] 
29] 
30] 
31] 
32] 
33] 
34] 
35] 
36] 
37] 



A. Suzuki et al., Nucl. Instrum. Methods Phys. Res. A 329 (1993) 299. 

F. Aharonian et al., Astronomy & Astrophysics 457 (2006) 899. 

J. Albert et al., Astrophys. J 674 (2008) 1037. 

F. Aharonian et al., Astrophys. J 614 (2004) 897. 

A. Borione et al., Astrophys. J 481 (1997) 313^ 



ht t p : / / umdgrb . umd . edu/ hawc / do cument s / H AWC- Sensi- 3 . p df 

F. Beacom & D. Kistler, Phys. Rev. D 75 (2007) 083001. 

S. Hunter et al, Astrophys. J 481 (1997) 205. 

S. LeBohec et al, Astrophys. J 539 (2000) 209. 

F. Aharonian et al., Astronomy & Astrophysics 375 (2001) 1008. 

F. Aharonian et al., Astropart. Phys. 17 (2002) 459. 

A. Borione et al., Astrophys. J 493 (1998) 175. 

M. Amenomori et al., Advances in Space Research 37 (2006) 1932. 

M. Amenomori et al., astro-ph/0611335vl. 



A. Abdo et al., Astrophys. J 688 (2008) 1078. 
T. Porter & R. Protheroe, J. Phys. G 23 (1997) 1765. 
N. Tateyama Sz J. Nishimura, in: Proceedings of the 28th International 
Cosmic Ray Conference, Tsukuba, 4, 2003, p. 2285. 
[39] V. Berezinsky et al., Astropart. Phys. 1 (1993) 281. 
[40] S. Wakely & D. Horan, http:/ /tevcat. uchicago.edu 



41] T. Totani & T. Takeuchi, Astrophys. J 570 (2002) 470. 



9 



Table 1 



Performance comparison between the Tibet AS+MD array with the CASA-MIA 
experiment (0) 





Energy threshold 


Angular 


\x det. 


AS array 


BG hadrons 




(75% trigger efficiency) 


resolution 


effective area 


effective area 


survival efficiency 


Tibet 




0.2° 




83000 m 2 


0.005% (100 TeV) 


~10 TeV 


10000 m 2 


AS+MD 




(100 TeV) 




(upgraded) 




CASA- 


~100 TeV 


2° 


2500 m 2 


230000 m 2 


2% (178 TeV) 


MIA 




(100 TeV) 






0.02% (646 TeV) 
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15 m 

Fig. 1. Schematic view of the Tibet AS+MD array. Open squares represent surface 
0.5 m 2 scintillation counters comprising the Tibet AS array. Filled squares represent 
52 m 2 pool cells of the MD array, set up beneath a 2.0 m thick soil and a 0.5 m thick 
concrete layer. The AS array shown here is an upgraded version of the present one, 
consisting of 1081 scintillation counters with an effective of 83000 m 2 . Enclosed 
by the dashed line is the fiducial area for selecting air shower events based on their 
core locations (See £j3]). 
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Fig. 2. Schematic view of a pool cell of the MD array (Left: top view, right: side 
view). The water volume is 7.2 m wide, 7.2 m long, and 1.5 m deep, with a 0.9 m 
thick air layer over it. Two downward- facing 20" 4> PMTs are installed on the ceiling. 
Two vertical pillars reinforce the structure. 
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Fig. 3. Quantum efficiency of 20" (j) PMTs assumed in the simulation. 
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Fig. 4. Distribution of the number of photoelectrons averaged over two PMTs (Npe) 
for vertically penetrating muons that uniformly hit the ceiling of an MD cell. 
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Fig. 5. Distribution of XiVpE as a function of Sp. T,p denotes the sum of particle den- 
sities for all the AS scintillation counters, and EiVpE the sum of iVpES from the MD 
pool cells. The green and grey dots correspond to the gamma- and cosmic-ray-in- 
duced air showers, respectively. The filled triangles (squares) and the accompanying 
error bars indicate values at the median, 20%, and 80% of the £iVp E distribution of 
gamma-ray (cosmic-ray) events in each Ep bin. The air showers producing no PEs 
are plotted at XiVpE = 0.5. The number of photoelectrons that each 20" (ft PMT 
can detect is assumed to saturate at ~300. The number of green and grey dots 
plotted do not reflect the true flux ratio of the gamma rays to the cosmic rays; 
for demonstration, a sufficient number of gamma-ray events are generated so that 
the high-energy gamma-ray events above 100 TeV can be recognized on the figure. 
Meanwhile, the fluxes of the gamma rays and the background cosmic rays under the 
assumption of the one- year observation of the Crab Nebula are considered in setting 
the SA^PE-based event-selection criterion, shown by the solid line, which maximizes 
the detection significance of the source. 
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Fig. 6. XiVpE distribution of gamma-ray and background cosmic-ray events in the 
10 TeV energy bin (Left; 100 < < 215) and in the 100 TeV energy bin (Right; 
1000 < Yl P < 2154), considering the fluxes of the gamma rays and the background 
cosmic rays under the assumption of the one-year observation of the Crab Nebula. 
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Fig. 7. The survival efficiencies of gamma rays (squares) and background cosmic 
rays (triangles) after the event selection by XiVpE as shown in Fig. [5j Note that 
the survival efficiencies defined here do not include the event reduction by the 
search window, which depends on the angular resolution. The background events 
are rejected down to the 2 x 1CP 3 level at 10 TeV, with 60% of the gamma-ray 
events remaining. The power-law fit of the cosmic-ray survival efficiencies (grey 
solid line) indicates that the background events are rejected down to the 5 x 10~ 5 
level at 100 TeV with 83% of the gamma-ray events remaining. The two points of 
background survival efficiency above 50 TeV were obtained from another data set 
of the background MC events generated above 10 TeV. 
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Fig. 8. T,p spectrum of the gamma-ray (open circles) and background cosmic-ray 
(filled squares) events expected, assuming the one-year observation of the Crab 
Nebula with a flux of F G = 6.45 x 1(T 10 x (£/0.3TeV)~ 2 - 6 [cnr" 2 s~ 1 TeV~ 1 ] and the 
integral flux of the cosmic rays of F CR (> 0.3TeV) = 1.03 x 10 -4 [cm~ 2 s —1 sr ]. The 
dashed (solid) lines show the spectrum before (after) the event selection by EiVpE 
as shown in Fig. [5j 
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Energy (TeV) 

Fig. 9. The 5a sensitivity of the Tibet AS+MD array to a point-like gamma-ray 
source for one-calendar- year observation (thick solid line), along with its 90% confi- 
dence level upper limit (thick dashed line). Note that the sensitivities above 100 TeV, 
where the number of background cosmic rays is fully suppressed down to less than 
one event, is defined as a flux corresponding to 10 gamma-ray events per year. 
For reference, the sensitivity of the present Tibet AS array is drawn by the thin 
dashed line. The Crab's flux points measured by HESS (0), MAGIC 0), HEGRA 
(j25l ). and Tibet AS (fl^ ) are represented by crosses, dots, open diamonds, and filled 
squares, respectively. The triangles with downward arrows represent 90% upper lim- 
its on the Crab's flux given by the CASA-MIA experiment (|26l ). For comparison, 
the 5<t sensitivities of currently operating H.E.S.S. as well as future experiments, 
HAWC and CTA, are plotted (0). 
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Fig. 10. The 5cr sensitivity of the Tibet AS+MD array to a spatially extended 
gamma-ray source for one-calendar-year observation (thick solid line). Note that 
the sensitivities above ~200 TeV, where the number of background cosmic rays is 
fully suppressed down to less than one event, is defined as a flux corresponding to 
10 gamma-ray events per year. The triangle represents the flux point of MGRO 
2019+37 at 20 TeV measured by MILAGRO A vr° decay gamma-ray flux 

calculated by Beacom and Kistler |2^ ) assuming the acceleration of protons up to 
1000 TeV is shown by the dashed line. 
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Fig. 11. The 5<r sensitivity of the Tibet AS+MD array towards diffuse gamma 
rays from the inner Galactic plane (left) and the outer Galactic plane (right) for 
a one-calendar-year observation (thick solid line). Note that the sensitivity above 
~1000 TeV is defined as a flux corresponding to 10 gamma-ray events per year. The 
EGRET data are represented by dots, assuming the differential spectral indices to be 
2.4(— 3.3) for the inner(outer) Galactic plane (0). Upper limits by Whipple (W) 



3C) with 99.9% C.L., HEGRA (H) (|3l|) with 99% C.L., HEGRA-AIROBICC (Ha) 
3 with 90% C.L., CASA-MIA (C-M) Q with 90% C.L., and the Tibet air shower 



experiment (34; 35|) with 99% C.L. are represented by downward arrows. MILAGRO 
(M) obtained a definite flux from the inner Galactic plane and an upper limit with 
95% C.L. from the outer Galactic plane at 15 TeV (1361 ). Theoretical curves of inverse 
Compton models by Porter and Protheroe (PP2.0 and PP2.4; dotted lines) (133), 
as well as those by Tateyama and Nishimura (TN2.0 and TN2.4; solid lines) (1381 ) . 
are included, where 2.0 and 2.4 are the source electron injection spectral indices 
they assumed. All the inverse Compton models assume an exponential cutoff of 
the source electron injection spectrum at 100 TeV. Theoretical curves arising from 
7T° — > 27 decay (BGHS) are given by Berezinsky et al (0). 
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Fig. 12. The energy spectra of 14 H.E.S.S. sources reported in QjJ) , along with the 5 
sigma sensitivity of the Tibet AS+MD array to a point-like gamma-ray source for 
one-calendar-year observation (black thick line) . The grey dots represent the integral 
fluxes calculated from the observed differential fluxes of the H.E.S.S. sources point 
by point, assuming their spectral indices. The grey thin lines show power-law fittings 
to the H.E.S.S. data points with exponential cutoffs assuming the cutoff energy of 
source protons to be 1000 TeV (28|). Note that the sensitivity above 100 TeV is 
defined as a flux corresponding to 10 gamma-ray events per year. 
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